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We present the Curci-Ferrari model on the lattice. In the massless case the topological interpre- 
tation of this model with its double BRST symmetry relates to the Neuberger 0/0 problem which 
we extend to include the ghost/anti-ghost symmetric formulation of the non-linear covariant Curci- 
Ferrari gauges on the lattice. The introduction of a Curci-Ferrari mass term, however, serves to 
regulate the 0/0 indeterminate form of physical observables observed by Neuberger. While such a 
mass m decontracts the double BRST/anti-BRST algebra, which is well-known to result in a loss 
of unitarity, observables can be meaningfully defined in the limit m — > via I'Hospital's rule. At 
finite m the topological nature of the partition function used as the gauge fixing device seems lost. 
We discuss the gauge parameter ^ and mass m dependence of the model and show how both cancel 
when m = m(^) is appropriately adjusted with ^. 

PACS numbers: 11.15.-q, 11. 15. Ha, ll.30.Ly, 12.38.-t, ll.30.Pb 
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I. INTRODUCTION 



In the covariant continuum formulation of gauge theo- 
ries, in terms of local field systems, one has to deal with 
the redundant degrees of freedom due to gauge invari- 
ance. Within the language of local quantum field theory, 
the machinery for that is based on the so-called Becchi- 
Rouet-Stora-Tyutin (BRST) symmetry which is a global 
symmetry and can be considered the quantum version 
of local gauge invariance [ll, y]- In short, one starts out 
from the representations of a BRST algebra on indefi- 
nite metric spaces with assuming the existence (and com- 
pleteness) of a nilpotent BRST charge Qb- The physi- 
cal Hilbert space can then be defined as the equivalence 
classes of BRST closed (which are annihilated by Qb) 
modulo exact states (which are BRST variations of oth- 
ers) . In QED this machinery reduces to the usual Gupta- 
Bleuler construction. For the generalization thereof, in 
non-Abelian gauge theories, all is well in perturbation 
theory also. Beyond perturbation theory, however, there 
is a problem with such a construction that has not been 
fully and comprehensively addressed as yet. It relates 
to the famous Gribov ambiguity p] , the existence of so- 
called Gribov copies that satisfy the Lorenz condition Q 
(or any other local gauge fixing condition) but are related 
by gauge transformations, and are thus physically equiv- 
alent. As a result of this ambiguity, the usual definitions 
of a BRST charge fail to be globally valid. 

A rigorous non-perturbative framework is provided by 
lattice gauge theory. Its strength and beauty derives from 
the fact that gauge- fixing is not required. However, in 
order to arrive at a non-perturbative definition of non- 
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Abelian gauge theories in the continuum, from a lattice 
formulation, we need to be able to perform the contin- 
uum limit in a formally watertight way. And there is the 
gap in our present understanding. The same problem 
as described above comes back to haunt us in another 
dress when attempting to fix a gauge via BRST formula- 
tions on the lattice. There it is known as the Neuberger 
problem which asserts that the expectation value of any 
gauge invariant (and thus physical) observable in a lat- 
tice BRST formulation will always be of the indefinite 
form 0/0 ;5,]. 

The BRST algebra requires the introduction of further 
unphysical degrees of freedom. These are the Faddeev- 
Popov ghosts and anti-ghosts which violate the Spin- 
Statistics Theorem of local quantum field theory on pos- 
itive definite metric (Hilbert) spaces. Contrary to what 
the name anti-ghost might suggest, however, in the usual 
linear covariant gauges the treatment of ghosts and anti- 
ghosts is completely asymmetric. On the other hand, it 
is also known for many years that it is possible to extend 
the BRST algebra to be entirely symmetric w.r.t. ghosts 
and anti-ghosts. This additional symmetry arises natu- 
rally in the Landau gauge but can also be extended to 
more general gauges, the so-called Curci-Ferrari gauges, 
at the expense of quartic ghost self-interactions. The 
most interesting feature of these gauges four our pur- 
pose, however, is that they allow the introduction of a 
mass term for gluons and ghosts 6] . While such a Curci- 
Ferrari mass m breaks the nilpotency of the BRST and 
anti-BRST charges, which is known to result in a loss of 
unitarity [T], |8[ and which therefore meant that this rela- 
tively old model received little attention for many years, 
it also serves to regulate the Neuberger zeroes in a lat- 
tice formulation. In [9] this was exemplified in a simple 
Abelian toy-model where the zeroes in the numerator and 
denominator of expectation values become proportional 
to m^ and allow to compute a finite value for rri^ — > 
via THospital's rule. 



For the SU{N) gauge theory on a finite four- 
dimensional lattice things are naturally much more com- 
plicated than in the toy model. In this paper we 
develop a full lattice formulation of the time-honored 
model by Curci and Ferrari with its decontracted double 
BRST/anti-BRST and ghost-mass term, as announced 
in [iO|- After introducing the general setup for double 
BRST on the lattice in Sec. II, we next review Neu- 
berger's no-go-theorem in a generalized version to include 
the ghost/anti-ghost symmetric case of the non-linear co- 
variant Curci-Ferrari gauges for m^ = in Sec. Ill, a 
case originally excluded by Neuberger. At non- vanishing 
Curci-Ferrari mass the partition function of the model 
used as the gauge-fixing device is shown to be polynomial 
in m'^ and to be thus non- vanishing, in a special gauge- 
parameter limit in Sec. IV. In this way regularizing the 
Neuberger zeroes, the leading power of that polynomial 
can be extracted from a suitable number of derivatives 
(w.r.t. rn^) before the limit m'^ — > is taken, in the spirit 
of rHospital's rule. This could provide a lattice BRST 
model without Neuberger problem. The massive Curci- 
Ferrari model is no-longer purely topological in nature, 
however, and as a result, its gauge-parameter ^ indepen- 
dence requires tuning of the Curci-Ferrari mass with ^ 
as explained in Sec. V. The gauge-orbit independence of 
this procedure is discussed in Sec. VI. A short summary 
is given in Sec. VII, and our conclusions and outlook are 
provided in Sec. VIII. 



II. DOUBLE BRST ON THE LATTICE 

For the topological lattice formulation of the double 
BRST symmetry of the ghost/anti-ghost symmetric co- 
variant gauges we start out from the standard gauge- 
fixing functional Vu[g] of covariant gauges which here 
assumes the role of a Morse potential on a gauge orbit, 

i^Retrt/,3,^. (1) 



^^[5]=-2^EEtrt/^- 



Here, in the first form, Uij £ SU{N) is the directed link 
variable connecting nearest neighbor sites i and j. The 
sum J ~ i denotes summation over all nearest neighbors 
J of site i. We assume periodic boundary conditions. 
The double sum thus runs twice over all links (ij), and 
with {// = Uji it is therefore equivalent to the simple 
sum over links in the second form, where U^.p, stands for 
the same link field U at position x in direction /i. The 
constant p is the normalization of the SU{N) genera- 
tors X. We use anti-Hermitian [X", X^] = f^bc^c .^^jth 
trX'^X'' — —pS""^. We exphcitly only need the funda- 
mental representation, where p = ptund = 1/2. 

As usual, under gauge transformations the link vari- 
ables U transform 

U.j ^Uf^^ glU,,g, . (2) 

BRST transformations s and anti-BRST transformations 
s in the topological setting do not act on the link vari- 



ables U directly, but on the gauge transformations gi like 
infinitesimal right translations in the gauge group with 
real ghost and anti-ghost Grassmann fields cf, c° as pa- 
rameters, respectively. 



sg = gX'^d' = gc, sg = gX'^if 



gc , 



(3) 



where we introduced Lie-algebra valued, anti-Hermitian 
ghost fields Ci = X'^cf with c] — —Ci, and analogous anti- 
ghost fields Ci = X^c". For consistency, we furthermore 
require 



sg 



t - (sn\'! - 



[sgr = -eg'' , sg^ = (sg) 



\t _ 



-eg 



(4) 



For the gauge-transformed link variables this then implies 

The BRST transformations for (anti)ghosts and Naka- 
nishi-Lautrup fields b are straightforward lattice ana- 
logues (per site) of their continuum counterparts, see, 
e.g., Refs. 0,111, 



sc- = -^{c X cr , 
s(f = 6'^-i(cxc)" 



3b^ = -i(cx6r~i((cxc)xc)" 



(6) 

(7) 
(8) 



The relatively obvious notation of using the "cross- 
product" herein refers to the structure constants for 
SU{N), for example, (c x cf = /°'"=c''c^ 

In the ghost/anti-ghost symmetric gauges as consid- 
ered here, the anti-BRST variations are obtained by sub- 
stituting c ^ c and c — > — c according to Faddeev-Popov 
conjugation. Thus, 



sc 

sc" 

sb" 



-b^-^{cxcr 
-^{cxcT, 



(9) 
(10) 



licxbr + l{{cx-c)xc)\ (11) 



The action of the topological lattice model for gauge fix- 
ing a la Faddeev-Popov with double BRST invariance can 
then be written in compact form as 



Sgf = i ss [Vu [g]+i—^tr 



(12) 



This is the lattice counterpart of the continuum gauge- 
fixing Lagrangian 

£^^ = ^ss{A''^A^^-i^ifc'') with 5gf = fd^xCov 

(13) 
in D Euclidean dimensions. 

For the purpose of a self-contained presentation we 
work out the double (anti-)BRST variation on the right 



of p^ explicitly in Appendix \X\ This leads to 

^GP - E { - 'b^Ft(U') ~^ctM,,t[c] (14) 

i 

+ ^bfbf + I (q X c,)2 } , 



where 



KiU') = -^ Etr(X'^(C/i^. -C/p) (15) 



2p 



defines, of course, the standard gauge-fixing form of co- 
variant gauges with the continuum limit, 



a^O 2 



Ff(C/9) "i:^ a'a^AS" + ©(a^) 



(16) 



The Faddeev-Popov operator Mpp" is obtained from the 
short-hand notation in ITJ" 



M, 



E< 






(17) 



and given explicitly for later reference in alternative 
forms in Eqs. (jAlip or (jAlSp . It is symmetric w.r.t. sim- 
lutaneous interchanges of color and site indices, and iden- 
tical to the one obtained in [iJl as the Hessian of Vu[g] 
from variations along one-parameter subgroups of the 
SU{N) gauge group. In the continuum limit it reduces 
to the symmetrized and thus Hcrmitian 



M, 



ab a^O 



vi(ai?«^ 



+ D^'d) S{x -y) + 0{a^) 



of the ghost/anti-ghost symmetric Curci-Ferrari gauges. 
In contrast, the Faddeev-Popov operator of the linear 
covariant gauges for ^ ^ is not a Hessian because it 
is not symmetric. It can be read off as a byproduct of 
our BRST derivation from Eq. (|A9[) . In particular, this 
non-symmetric Faddeev-Popov operator needs to be used 
when implementing other linear-covariant gauges such as 
the Feynman gauge with ^ == 1 on the lattice as discussed 
in [M, T5] . In Landau gauge C = the distinction is an 
illusion. To keep the symmetric Hessian for ^ 7^ 0, how- 
ever, is only possible within the ghost/anti-ghost sym- 
metric framework where it necessarily comes along with 
the quartic ghost self- interactions in p4)) . 

The full symmetry of the ghost/anti-ghost symmetric 
Curci-Ferrari gauges [g, Il2l | is given by a semidirect prod- 
uct of a global S'L(2,IR), which includes ghost number 
and Faddeev-Popov conjugation, with the BRST/anti- 
BRST symmetries as used above ^]. This is the global 
symmetry of the Landau gauge, and it is sometimes re- 
ferred to as extended BRST symmetry, see [l|. 

Among the general class of all covariant gauges pTJ . 
with a Lagrangian which is polynomial in the fields, 
Lorentz, globally gauge and BRST invariant, and renor- 
malizable in D = 4, the ghost/anti-ghost symmetric case 
is special and interesting in that it allows to smoothly 
connect to the Landau gauge for ^ — > 0, without chang- 
ing the global symmetry properties. 



In particular, introducing with [llj a second gauge pa- 
rameter P G [0, 1], to interpolate between the various 
generalized covariant gauges, the linear covariant gauges 
of standard Faddeev-Popov theory correspond to the line 
/3 = in the two gauge-parameter plane (^, (3). Along 
this line, the global symmetry changes abruptly when 
reaching the Landau gauge limit; and for P — 1, one 
obtains a mirror image of standard Faddeev-Popov the- 
ory with the roles of ghosts and anti-ghosts interchanged. 
The ghost/anti-ghost symmetric gauges discussed here 
then correspond to the line P = 1/2. The £, — gauge is 
/3-independent. The whole interval for /3 g [0, 1] at ^ = 
is equivalent and corresponds to the Landau gauge. The 
important difference is, however, that the SL{2, IR) sym- 
metric line at /3 = 1/2 provides a unique class of covariant 
gauges which share the full extended BRST symmetry of 
the Landau gauge for any value of ^. The limit ^ — *■ is 
thus a smooth one, as far as this symmetry is concerned, 
only along the line of /3 = 1/2. The price to pay are the 
quartic ghost self-interactions in (I14p which again vanish 
only in the Landau gauge limit. 

For a further discussion of the general ghost creat- 



ing gauges, and their geometrical interpretation, see 12 1 
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The one-loop renormalization was first discussed in 
for explicit calculations of renormalization constants and 
anomalous dimensions of the ghost/anti-ghost symmetric 
case up to including the three- loop level, see @, [l^| . The 
Dyson-Schwinger equations of these gauges were studied 
in [l7|. A non-renormalization theorem relating to the 
Curci-Ferrari mass was recently reported in [l8| . 



III. THE NEUBERGER PROBLEM 

Following Neuberger, we introduce an auxiliary param- 
eter t in the Euclidean partition function to be used as 
the gauge-fixing device via the Faddeev-Popov procedure 
of inserting unity into the unfixed partition function of 
SU{N) lattice gauge theory. The gauge- fixing action of 
the double BRST invariant model given by p^ consists 
of two terms both of which are separately BRST (and 
anti-BRST) exact. Multiplying the T"' term in ^ by 
the real parameter t amounts to a mere redefinition of 
the Morse potential which should have no further effect. 
We can therefore write the gauge-fixing partition func- 
tion with double BRST, 



ZaAt) = / d[g,b,c,c] X (18) 

exp I - issyt Vu[g] + *^ E ^'' ^'^V } ' 

i 

which is independent of the set of link variables {U} and 
the gauge parameter ^ because of its topological nature. 
Moreover, the t independence is really no different from 
the ^ independence here, and it is thus rather obvious. 
Explicitly, the derivative with respect to t (or ^) produces 
the expectation value of a BRST exact operator which 



vanishes, 



z'oAt) = 0. 



(19) 



At i = on the other hand, we obtain with the BRST 
variations given in (jA4p and (jA12[) of Appendix \^ 

Zgf(O) = JV f d[b,c,c] X (20) 

i 

where the volume of the gauge group on the lattice, from 
the invariant integrations Yii dgi via the Haar measure 
over gi g SU{N) per site i, is absorbed in the con- 
stant Af. The Gaussian integrations over the Nakanishi- 
Lautrup fields b are also well-defined and produce a factor 
{2T:/^)(^"-^y^ per site. 

One might be tempted to conclude at this point that 
the quartic ghost self-interactions in (P0|) might remove 
the uncompensated Grassmann integrations of the linear 
covariant gauges where no such self-interactions occur. 
The ghost/anti-ghost integrations at i = also factor- 
ize into independent integrations dc^dc^ over 2(A^^ — 1) 
Grassmann variables per site. For TV = 3, for exam- 
ple, the 4**^ order term of the exponential in ([^0]) pro- 
duces a monomial in c° , cf which contains each of these 
16 Grassmann variables exactly once, so that their in- 
tegration might produce a non- vanishing result. This is 
not the case, however. Working out the prefactor of this 
monomial, as we will do explicitly in the more general 
case with including a non-vanishing Curci-Ferrari mass 
m below, one finds that the prefactor of this term in ([^D]) 
vanishes in the massless case and thus. 



topological invariant, has of course no effect in terms of 
fixing a gauge. In the present formulation, with Zgf(O) 
in (|20p , the independent Grassmann integrations per site 
of the quartic-ghost term which contains the curvature 
of SU{N) each compute its Euler characteristic via the 
Gauss-Bonnet theorem [2J| . This explicitly produces one 
factor of zero per site on the lattice. And it provides the 
topological explanation for the vanishing of the prefac- 
tor of the corresponding monomial of degree 2(A^^ — 1) 
in the Grassmann variables c, c, which could otherwise 
exist in the expansion of the exponential in (|20p for all 
odd N . For A'^ = 3, for example, the zero in this pref- 
actor arises, upon normalordering, from a cancellation of 
368 non- vanishing individual terms when expanding the 
square of the square of the quartic ghost self-interaction. 
This cancellation would be rather unnatural to arise ac- 
cidentally, without such explanation. 

The vanishing of the gauge-fixing partition function at 
t = part in Neuberger's argument, in the ghost/anti- 
ghost symmetric gauges with their iS'L(2,IR)x double 
BRST symmetry, therefore most directly refiects the 
topological origin of the Neuberger zero. Eq. (f20|) pre- 
cisely represents a product of one Gauss-Bonnet integral 
expression for x{SU{N)) per site of the lattice. 

Note that the gauge parameter ^ can be removed com- 
pletely from the expression for Zgf(O) in Eq. (|20p by a 
rescaling v^ 5 — > 6 and -^ c^ c, ^^ c ^ c, which leaves 
the integration measure unchanged. The same rescaling 
for the full gauge- fixing partition function Zo^{t) in (fTSj) . 
which amounts to replacing the action in S'gf in (|14p by 

So.{t) = E { - 'i^1P^iU') - zic?MFP?[c] (22) 



^gf(O) - 



(21) 



:fe?6? 
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(Cj X C^f I , 



Because of the t- independence (|T9)) , this implies the van- 
ishing of the gauge-fixing partition function P^ of the 
ghost/anti-ghost or SL{2, IR) symmetric formulation with 
double BRST invariance in the same way as that of stan- 
dard Faddeev-Popov theory observed in Q. As for the 
latter, the sign-weighted sum over all Gribov copies, as 
originally proposed to generalize the Faddeev-Popov pro- 
cedure in presence of Gribov copies jlE, lSA], vanishes. 

This cancellation of Gribov copies is well-known [2ll |. 
The fact that it also arises here, in the ghost/anti-ghost 
symmetric formulation with its quartic self-interactions, 
directly relates to the topological interpretation [22, [23| 
of the Neuberger zero: Z^-p can be viewed as the partition 
function of a Witten-type topological model to compute 
the Euler characteristic x of the gauge group. On the 
lattice the gauge group is a direct product of SU{Nys 
per site, and because the Euler characteristic factorises. 



Zn 



X(SU{N)*''''=') = x(S'C/(A^))#''*"' = 



:^sites 



For t — Q the action in p^ decouples from the link-field 
configuration and Zgf(O), albeit computing the same 



furthermore shows that t and ^ really represent a single 
parameter i/VJ- Setting t = in Neuberger's argument 
is therefore the same as the ^ — > oo limit which is usually 
what is considered as the Gauss-Bonnet limit in topologi- 
cal quantum field theory [2J]. As mentioned above, there 
is no gauge-fixing in this limit, but it provides a simple 
way to compute the value (zero here) of the partition 
function which is independent of t/^/^. 

In the opposite limit, that of the Landau gauge ■f ^ 
01' W? ^ oOj of course, Zoj,{t) still reduces to the sign- 
weighted sum over all Gribov copies as usual [1^, l20| , 



ZoAt) 



J2 sign(detAfFp(C/9*")) 



(23) 



{3(''} 



which because of the t (and £,) inde pen dence ( T9|) thus 
computes the same topological zero [2ll . |22| . l23j , in this 
case via the Poincare-Hopf theorem [2J| . 



IV. THE MASSIVE CURCI-FERRARI MODEL 
ON THE LATTICE 

In the previous section we have seen that the quar- 
tic ghost self- interactions of the SL{2, IR)xi double BRST 
symmetric Curci- Ferrari gauges have no effect on the dis- 
astrous conclusion of the 0/0 problem in lattice BRST. 
They rather serve to reveal most clearly the topological 
origin of this problem. 

We will demonstrate explicitly below that this zero can 
be regularized, however, by introducing a Curci-Ferrari 
mass m, as proposed in [^, [I3|- The gauge- fixing action 
5'gf is thereby once more replaced by 



tVuig] + ^e^trQc) (24) 



S'„.gf(0 = «(«■»" 



(where we dropped in the 2""^ term the factor l/(2p) = 1, 
in the fundamental representation). The BSRT and anti- 
BRST transformations of C/^, c and c in Eqs. (O, (O, 
^ and (O, (Uni) of Sect. |TT] remain unchanged. Those 
for the Nakanishi-Lautrup 6- fields, Eqs. ^ and (|lip . are 
replaced by |13 |. 

s6" = im^c" -i(cx6)"-i((cxc)xc)", (25) 

sb" = im^ c" - Uc X hf + i((c x c) x c)" . (26) 
2 8 

In the derivation of the explicit form for S'„,GF(i), us- 
ing these modified (anti-)BRST transformations, the only 
modification in comparison to Sect. HIl and App.EI arises 
from s(c°6") in (jA12[) . which now becomes, 

s{eb'') = -im^ (fe + b^b" + i{cxcf . (27) 

The additional first term on the right contributes an 
additional term — i(^/2)TO^c"c" to the gauge-fixing La- 
grangian, c.f., Eq. (jASp . Together with the same contri- 
bution from the explicit mass term —i{^/2)m'^cfcf in (IM)) 
we therefore obtain twice that as the final ghost mass- 
term of the massive Curci-Ferrari model (this subtlety 
will be worth remembering for later) . The action of the 
massive Curci-Ferrari model therefore becomes, explic- 

itly, 



S^cAt) = mHVu[g] + 



(28) 



E 



{ - itb°;F^{U3) ~ itc1M^^1[c 



+ 1 6? 6? - im 



tm^^c^cf + |(c. xc.)^} 



BRST and anti-BRST transformations are no-longer 
nilpotent at finite m^ , but we have [l|, @, [l2] 



2 ■ 2 + 

s ~ tm a 



-2 - 2 - 

s = —tm a 



ss + ss = —im a 



(29) 



where a^ and a^ generate the global SL{2, IR) in- 
cluding ghost number and Faddeev-Popov conjugation. 



The Curci-Ferrari mass decontracts the sl{2, IR)x double 
BRST algebra of the massless case to the osp{l\2) super- 
algebra extension of the Lie algebra of the 3-dimensional 
Lorentz group SL{2, IR). Conversely, the m^ -^ limit is 
interpreted as an Inonu-Wigner contraction of the simple 
superalgebraosp(l|2) [l|,[i2|. The BRST and anti-BRST 
invariance of the massive Curci-Ferrari action in (|24p it- 
self follows readily from this algebra as given in (|29p . 
noting that only c and c transform non-trivially under 
theS'L(2,IR). 

We emphasize that this algebra decontraction has from 
the very beginning been known to lead to a breakdown 
of unitarity when attempting a BRST cohomology con- 
struction of a physical Hilbert space in analogy to the 
massless case [6|. In fact, explicit examples exist for 
states of negative norm surviving in any such construc- 
tion [7], Q . They do not belong to BRST quartets and 
can therefore not be removed by the quartet mechanism 
l]. Only through the algebra contraction by m^ — > do 
these states reduce to zero norm components which have 
no effect on the physical S'-matrix elements. 

Here we deliberately do not want to interpret the mass 
parameter by Curci and Ferrari as a physical mass. It 
rather serves to meaningfully define a limit m^ — > on 
the lattice, perhaps in parallel with the continuum limit, 
to recover nilpotent (anti-)BRST transformations. 

To study the parameter dependence, we first define 
the partition function of the massive Curci-Ferrari model, 
explicitly listing all three parameters (even though these 
again really only represent 2 independent ones as we will 
show below), 

^.„GF(t,^,m2) ^ J d[g,b,c,c] exp{ - S„,aAt)} , (30) 

with S'„,GF(i) from P^ or (|^ . We note in passing that 
the terms proportional to m? in the massive Curci-Ferrari 
action ((28|) are given by 



o(t,e)^iPc;[.9]-*eEcX , 



(31) 



or, in the continuum, 

O{t,0=Jd''x[^Alix)A'^^ix)~t^cAx)c'^{x)) . (32) 

For i = 1 this coincides with the on-shell BRST invari- 
ant (at m^ — 0) operator proposed by Kondo as a pos- 
sible candidate for a dimension 2 condensate [2^. The 
doubling of the explicit ghost mass-term in ([12]), by the 
BRST variation of cb in (P7|) as mentioned above, is cru- 
cial here. Without this difference in the relative factor of 
2 between the two terms in ©(t, ^) and the gauge fixing 
functional 



tW^,^tVu[g]~i^Y.^^''^ 



(33) 



one could not have both, the on-shell BRST invariance 
of O and the gauge-fixing action in p2|) from the double 
BRST variation S'gf = ss Wgf j at the same time. 



The observation that the mass terms in (pS)) are given 
by m?0[t, ^) could in principle be used to obtain the ex- 
pectation value of Kondo's operator from the derivative 



{o{t,i))^ 



d 



dr 



lnZ„,GF(i,C;"i^ 



(34) 



upon insertion into the unfixed partition function of lat- 
tice gauge theory, i.e., with taking the additional expec- 
tation value in the gauge-field ensemble. As any other 
observable at m^ = this expectation value as it stands, 
unfortunately, of course also suffers from Neuberger's 0/0 
problem of lattice BRST. 

In order to demonstrate that the Curci-Ferrari mass 
regulates the Neuberger zero, for t = we will verify by 
explicit calculation that 



Z„oF(0,e,™') ^ 
In fact, from ((30l) . (IM]) . 

Z,„GF{0,^,m^) = Af d[b,c,c] x 



(35) 



(36) 



exp { - ^ (i6^6^ - tm'^c'^c'^ + f (c. x c,)') } , 

i 

which again factorises into independent Grassmann (and 
6-field) integrations per site on the lattice. Using the 
same rescaling y^b -^ b and -\/Jc — > c, y^c — * c as 
mentioned in the last section, we obtain. 



(y^(2^)(^^'i)/2/^(m2v^) 
where Vn is the group volume of SU{N), and 

lN{fn^) = / W d[i^)dc'' 



(37) 



9^sites 



(38) 



'{^ 



exp < im c-c 
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(ex c 



'1 



where we used the rather obvious abbreviations c- c — 
c^c", (c X c)° = f'^'c^c", and m^ = m'^^/l. Note that we 
define the Grassmann integration measure to include the 
imaginary unit i with the real anti-ghosts c so as to re- 
produce the result of integrating over complex conjugate 
Grassmann variables c°'±ic°' . Expanding the exponential 
and collecting the relevant powers in the ghost/anti-ghost 
variables, for SU{2) we straight-forwardly obtain. 



hifh^) 



('+H 



3-2 /-, , 4_4> 

-TO 1 H — m 
4 V 3 



(39) 



For SU{3) the computation is a bit more tedious, the 
result is 



h{m ) = —to'' 1 + 4™" -I- —TO** -I- —to' 
64 V 15 45 



(40) 



In both cases we factorised the leading power for to^ -^ 0. 
lN{fn^) is polynomial in to^ = vn? -^/J of degree N"^— 1 , for 



all N. The successively lower powers of fh? decrease by 
2 in each step in this polynomial, reflecting an increasing 
power of the quartic ghost self-interactions contributing 
to each term. Therefore, the polynomials /Ar(TO^) are 
odd/even in fh? for N even/odd. 

Because the polynomial is odd for all even N there 
can thus not be an order-zero term in the first place. 
The powers of the quartic interactions alone never match 
the number of independent Grassmann variables, and the 
Neuberger zero at to^ = arises rather trivially for even 
iV, for the same reason that the Euler characteristic of 
an odd-dimensional manifold, here of dimension A^^ — 1, 
necessarily vanishes. 

For N odd, Iiq{'fh?) is an even polynomial which could 
in principle have an order zero, constant term. The fact 
that this term is absent, e.g., as explicitly verified for 
5J7(3) in (|40| . reflects the vanishing of the Euler char- 
acteristic of SU{N) also for odd N , as mentioned above. 
The in this case even dimension N"^ — 1 of the algebra is 
irrelevant here, because, for the purpose of cohomology, 
the parameter space of SU{N) behaves as a product of 
odd-dimensional spheres ^^ x 5^ x 5^ x • • • S'^^''^ 0. 

The polynomials lN(jn^) do not have a constant term 
in either case and therefore vanish with ffi^ —^ 0, i.e., 
-fAf(O) = 0, as expected. Moreover, the scaling argument 
used here and in the last section shows that the partition 
function (|30p of the massive Curci-Ferrari model can only 
depend on two of the three parameters. 



.(i,C,m2) = f{t/^,^m^) 



(41) 



An independent route of deriving this generic form, from 
the equations of motions, will be presented below. In 
this section we explicitly obtained /(O, y) with y — fn^ 
to constrain this function f{x,y) of two variables along 
the X ~ t/y/^ = line, and verified that 



^„GF(0,e,m2) = /(0,em4)cx 



(^^4)#sitcs/2^ AT =2 



for m^ -^ 0. Because of the topological explanation of the 
zero obtained in this limit, i.e., /(0,0) = 0, as discussed 
in the last section, this actually constrains / to vanish 
along the entire y = line, f{x, 0) = for all x = tj \fi,. 
For X = we can in principle therefore define a non- 
vanishing, finite limit. 



lim (^TO, ) 



-Wtot 



j(0,<^,TO^) — const. 



(42) 



with an appropriate power A^tot = =^ of sites on a finite 
lattice for odd A^, or half that for even N . This constant 
could thus be inserted into the unfixed lattice gauge the- 
ory measure without harm, i.e., avoiding the zero in (j21|l . 
Because x = t/y^ = 0, however, this still has no effect in 
terms of gauge- fixing by the Faddeev-Popov procedure ei- 
ther. We need to get away from a; = 0, at least by a small 
amount, in order to suppress those parts of the gauge or- 
bits with large violations of the Lorcnz condition. At a 



finite Curci- Ferrari mass m?, however, this is aggravated 
by the fact that the gauge-fixing partition function of 
the Curci- Ferrari model is no-longer that of a topological 
model, and we thus no-longer have the i-independence 
(or x-independence) of P^ either. We can therefore not 
as yet conclude at this point that the constant in (|42|) will 
essentially remain unchanged when going to some finite 
X = t/VS. 7^ as we must. 

We are not quite there yet, and we will therefore have 
to have a closer look at the parameter dependence of the 
massive Curci-Ferrari model in the next section. 



V. PARAMETER DEPENDENCES 



From Eqs. ^^ and ([M]) or (gS]) we immediately obtain 
the following (logarithmic) derivatives, 

d 
d 



2e^lnZ„oF(i,e,"''') 



i{ [ss — im' 



'■){-^^Y.^4))r 



m 



dm'-' 



lnZ„oF(i,C,m2) = -{m^Oit,0), 



(43) 



where the subscripts m^ on the right denote expectation 
values within the Curci-Ferrari model at finite mass. In 
particular, the derivative w.r.t. m^ in the last line differs 
from ([M]) only in that m^ has not been set to zero here 
yet. All these expectation values can, in general, depend 
on the link- field configuration {[/} which acts as a back- 
ground field to the model. Independence of {U} is only 
guaranteed to hold in the topological limit ni^ -^ 0. 
From the definition of O in (I^TI) , we thus find that 



(44) 



4 + 24-m^^)lnZ_(^,C,m^) 



-z(..-O(t,0),„. . 

The standard argument that the expectation value of an 
(anti-)BRST exact operator vanishes does not hold at fi- 
nite m^. Neither are BRST and anti-BRST variations 
nilpotent, nor is O invariant under the BRST or anti- 
BRST transformations. However, the equations of mo- 
tion for (anti-)ghost and Nakanishi-Lautrup fields on the 
lattice, i.e., their lattice Dyson-Schwinger equations, can 
be used to show that, indeed. 



(..-o(t,o),„.-o, 

even at finite m^. Therefore, 



(45) 



, d , r,y- d 2 



Z„^aAt,tm^) = . (46) 



This differential equation entails that we can write the 
partition function of the model in the generic form (|¥l]l . 



As we already did in the previous sections, we therefore 



continue to use the new parameters x = t/\/^ and 



^V^ from now on, writing 



,(a;,m^) 



(47) 



Again using rescaled fields y^b — > b, -^c — > c, v^c ^ c 
and with </Js -+ s, ^s -> s, so that the (anti-)BRST 
transformations of Eq. Q ~ (jlip remain formally un- 
changed, the only modification is the replacement of m^ 
by m^ in those of the massive model in Eqs. ((25l) . (|26p . 
Correspondingly, all other relations above are then con- 
verted by the formal replacements S, —> I, t —>■ x and 
m^ —> fn? . In particular. 



'-'iiiGfI.^;) 



^{ss-^fh')(xVu[g]-'-Y,c^c^ (48) 

i 
i 



with 



0{x) =xVu[g] ~iJ2^^ 



(49) 



The two independent derivatives left, are readily read off 
in an analogous way to give 



-- lnZ,„GF(a;,TO^) 
ox 

liiZ^Gp{x,m^) 



dfn? 



'i{{ss~im^)Vu[g\)- 



(50) 



In absence of a topological argument for the gauge pa- 
rameter independence at finite Curci-Ferrari mass, the 
best we can do to achieve independence of a; = VV? is 
to allow an x dependent mass parameter m? = fn?{x). 
In particular, the a; = results of the previous section 
are then to be interpreted as being expressed in terms of 
TO^(O). These results will remain unchanged for a; ^ 0, if 
we adjust the mass function fh?{x) with x in the partition 
function Z,„gf, accordingly That is, if 







-Z„oF(x,m 



{x)) 



(51) 



d 



dm d 



= \7^ + ^^l^] Z.^o.{x,m\x)) 



dx 



dx dm 



From Eqs. ([50)1 we see that this requires that 



dfn? 
dx 



X{ss-im^)Vu[g]), 



(52) 



This might not appear to be a very profound insight, 
because we simply arranged matters by hand to achieve 
gauge-parameter independence in this way. The crucial 
question at this point is, whether the tuning of the Curci- 
Ferrari mass parameter with x is possible independent of 
the link configuration {[/} which is far from obvious here. 
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Otherwise we would have to choose a different trajectory 
in the parameter space (a;,m^) for different gauge orbits 
which would be of little use then, as far as the Faddeev- 
Popov gauge-fixing procedure is concerned. If it is possi- 
ble, on the other hand, we can then use the value of the 
mass ttiq = m^(0) at x = to regulate the Neuberger 
zero and use the x and {f/} independent, non-vanishing 
and finite constant 



lim {mty^'°'Z^ 



t{x,m^{x)) ~ const. 



(53) 



as the starting definition of Faddeev-Popov gauge fixing 
on the lattice. Then, of course, we would also expect that 
there should be a topological meaning to this constant 
which is so far, however, unfortunately unknown to us. 



VI. ORBIT INDEPENDENT GAUGE- 
PARAMETER EXPANSION OF THE CF MASS 

As we have seen in the previous section, the gauge- 
parameter independence of the gauge-fixing partition 
function Z,„gf of the massive Curci- Ferrari model will in 
general require the rescaled Curci-Ferrari mass parame- 
ter fh? — m}\J\ to depend on the gauge parameter in a 
non-trivial way, i.e., ff? = m^(x), via x — t/^/^. 

Gauge-parameter independence requires the derivative 
of fh?{x) to be given by Eq. ([52|) . Together with the 
condition that m^(0) = TOq we can use this equation to 
obtain the coefficients of iv? (x) in a Taylor series expan- 
sion around x = 0, where we can do explicit calculations. 
Importantly, we can then verify that these coefficients 
will not depend on the gauge orbit, i.e., on (the class of 
gauge-equivalent) link configurations {U}. Being based 
on the tensor method of invariant integrations over the 
gauge group elements per site i of the lattice, this pertur- 
bative expansion of the Curci-Ferrari mass parameter at 
small X will in fact be gauge-orbit independent at every 
order. As always, of course, nothing can be learned in 
such an expansion about possible non-analytic contribu- 
tions. We therefore assume the analyticity of the mas- 
sive Curci-Ferrari model in the 'would-be-Gauss-Bonnet' 
limit a; — > 0. This should surely be valid in the massless 
limit, but we need to assume here that the limits a; — > 
and TOq — > can be interchanged, in addition. While this 
is all well on a finite lattice, it certainly needs to be kept 
in mind when studying the model in the infinite volume 
and continuum limits. 

On a finite lattice, it is relatively straightforward 
to show that w? is in fact independent of the gauge- 
parameter X at l*** oder in this expansion, i.e., that 



dx 



x=0 



(54) 



This is simply because the numerator in Eq. ((521) vanishes 
at a; = while the denominator is a finite number. To see 
this explicitly, first consider at a; = Eq. ([36]) with our 



new variables and rescaled fields, before the gauge-group 
integrations. 



Z,„gf(0,TOo) =: / d[g,b,c,c] x 



(55) 



ITTln C< 



■i' Ci ] r 



As mentioned above, it decouples from the link-field con- 
figuration and factorises. Relative to this partition func- 
tion, we obtain. 



:<c5 



N^-l In {ml) ' 



(56) 



which is easily verified from the rules of Grassmann inte- 
gration and Eqs. (I36])-(I381). Because 0(0) = -i^^i cfcf, 
the denominator in (j52p at a: = is obtained from the 
trace in (l56l). 



(0(-)), 



= {-^Y.- 



— (#sites) X 






In [ml) ' 



which on a finite lattice with non-vanishing mp is finite. 
For the numerator in (|52p we have to compute 



— lnZ„,GF(a;,"i^) = -*( (i 



')Vu[9]), 



{i{b, F) + i{c, Afppc) - m^Vu[g] ), 



(58) 



where the brackets are introduced for summation over 
site and color indices. At a; = 0, the 6-field integration 
is Gaussian and the first term in ((58l) . linear in 6, there- 
fore vanishes. Because the gauge fields decouple from 
the measure in ([55)) , Eq. (|56p produces the trace of the 
Faddeev-Popov matrix Afpp in the second term on the 
right of dSgi for a; == 0. With 



trAfpp([/9 



-2C'^Vu[g] 



(59) 



where Cj is the quadratic Casimir invariant in represen- 
tation p, 6*2 = {N—l/N)/2 in the fundamental one, we 
see that the second becomes proportional to the third, 
and both proportional to the expectation value of the 
Morse potential Vjj [g] . This is linear in U^ and contains 
exactly one element g^ (or gj) in each and every term for 
which the invariant integration over the gauge group at 
the particular site i produces a zero. Therefore, 



{Vu[g]), 



x=0 



0, 



(60) 



which establishes (|54l) . It means that the derivative of 
fh{x) is of the order x near a; = 0, or 



m2(a;) = ml + Oix^) 



(61) 



In order to compute the constant in the second order 
term, and verify that it is non- vanishing and independent 



of {U}, we can consider the derivative w.r.t. x of the 
numerator in (1521). or 






InZ, 



[x^w?) = 
((z(6,F)+i(c,Mp 



(62) 



m^Vr 



[5])'), 



- ( i(6, F) + i(c, Afppc) - mVa[g] ) -, . 

The second (disconnected) term again vanishes at 2; = 0. 
Expanding the square in the first, we once more use that 
every term hnear in the 6- field will also vanish at x = 0, 
and therefore, 



& 



^lnZ„op(x,m2) ^ = {{^b,F) ) 



dx 



'm2 



+ ((lC,Mppc)^)^, 



(2,0 



ff^l{vS[9])f 



(63) 



l2. 



-2fhl{{ic,M^^c)Vu[g]) 



m^ ' 

where we used shorthand (...), 7^2 q = (... )m2 U=o- 

We calculate and discuss each of the individual terms 
on the right of (p5|) separately in Appendix [B] (with the 
exception that, by the same argument that led to ([55)l . 
the last two are essentially the same again, i.e., both 
ex ( V^ )m2 g). In particular, we show there that they are 
indeed all independent of {U}. The results for SU{N) 
gauge-groups in D Euclidean dimensions arc summarized 
in Table [H where Im is defined in Eq. p8)l and Jn anal- 
ogously by 



JN{rri 



f lldiz^dc^ [-licxcf) (64) 

-' a=l 

I < ifh c-c — -{c X c)^ > . 



exp < im c-c 

A comparison of the integral expression for J^ with the 
analogous one for /jv in Eq. psp shows that with an 
explicit interaction inserted in the integral, Jisjirri?) is a 
polynomial in fh? of two orders less than Ipf(fn^). Just 
as /jv(^^), this polynomial has no zero-order, constant 
term because the Euler characteristic of SU{N) vanishes 
and we essentially obtain a Gauss-Bonnet integral for 
fh^ — + again. Therefore, 

Jjv(to^) ~ /Ar(TO^) ^ , for m^ 
Exphcitly, for SU{2), see Appendix [Bl 



(65) 



J2{fh') = ^l!^{fh') 



3 ^2 
- m 
4 



(66) 



while for SU{3) a straightforward but a bit more tedious 
computation analogous to that used for Eq. ([^0]) yields. 






3 m' 



(67) 



With the results in Table |T] we can go back to the 
derivative of the Curci- Ferrari mass parameter fh?{x) 



w.r.t. X. Recall that the independence of the massive 
Curci-Ferrari model on the gauge parameter x = i/V^ 
requires this derivative to be given by the ratio, c./., 
Eqs. dnH) and (E 



dx 



d_ 
dx 



lnZ„GF(a;,"i^ 



{0{x)), 



(68) 



With ([63|) and the results from the table, in the limit 
X — > we therefore find, 



> x2D < 

dx L 



N^ 



iN(nio) 



N r^{ml) 



for A^ > 3, and 

dfh x^o 






2m 



(69) 
Oivl • 



^2i?(3i» 

4 Hml) 



dx L" mfhl) 4 I^ifhl) 

'ojTpiy - 2mg} = -x2Dmlj2{ml) 

for S'C/(2), where, using Eq. ^ for /a, 
3 + 18ml, + 8ml, 



(70) 



72(^0) = 



3-|-12m(; 



(71) 



In either case, the expansion of the Curci-Ferrari mass 
parameter around a; = for SU{N) in D dimensions can 
finally be written in the general form 



m (x) 



x^DmljNifhl) + 0{x^) , (72) 



where jn{ttiq) is a ratio of polynomials in ttiq obtained 
via dMl), for A^ > 3, from 



iNiml) = fF 1 + 



N 






N^ - 1 



(73) 



N -'V mg/^(mg) 

For A^ = 3, using (PD|) and (|67p . we find explicitly. 



73 ("Iq) 



1 



(74) 



90 + 855 mi 



1692 m§ + 1088 mj^ _^ 192 ^J^ 



90 + 720 m^ -h 1152 a« + 512 aj^ 

To work out 7jv for general A^ in the limit TOq -^ 0, 
remember that lN{m'^) is an odd/even polynomial in to§ 
for A^ even/odd. In either case there is no constant term, 
/7v(?tiq) -^ for mf, -^ 0. Therefore 



In {ml) ffil^o 



mVNi^l) 



1 , A^ even, 
i , A^ odd. 



(75) 



For the leading power n of fn^ in Ii^{m1) near m§ = 0, 
which is n = 1 when N is even and n = 2 when A^ is 
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Expectation Value 



N = 2 



iV> 3 



<(ic,Mppc)\/^[3])^,,^ 



-2D (#sites) 3 
2D (#sites) 2 

3 + 12m(S 



-2D (#sites) 

-2D (#sites) 

— (#sites) 



3mg -I- 4m[; 



6mo 



3mg -I- 4mg 
3 + 12m(5 



3mi5 + 4mg 



-2D(#sites) iV 



2D (#sites) 

1 i'N{ff>:i) 



-2D (#sites) 

-2D (#sites) 

— (#sites) 



N In (ml) 
2_ JNJml) 
N In (ml) 
I'N(ml) 
In (ml) 



TABLE I: The individual terms in Eq. (|63p for the numerator of the derivative in (|52|l in 
the limit a; — > together with the leading contribution to the denominator from Eq. (|57|l . 



odd, we need to expand the exponential of the ghost self- 
interactions in its integral representation psp to a power 
p such that 

n + 2p = iV^ - 1 . 

Otherwise the Grassmann integrations over A^^ — 1 ghost 
(and anti-ghost) variables will produce zero. Therefore, 



N^ -2 



N^-3 



N even, 



TV odd. 



(76) 



Comparing the definition of Jn in Eq. (jB37|) to that of 
In in ([55]) , we see that the exponential of the ghost self- 
interactions in the integral representation of Jn needs to 
be expanded to one power less, i.e., to the power p — 1 
for the leading term. Comparing the prefactors of these 
terms in each case we therefore find. 



JN(fhl) m^-O p! 



iNifhl) (p-1)! 

With dTHl) and ^^ this then imphes. 



(77) 



N^ 



mpNiml) 



2 
TV^ -3 



N even, 
N odd. 



(78) 



For even A^ > 4 we thus obtain from Eq. ([731 



7w 



1'^ 



iV^-2 



N^ 



N 



1 

iV ' 



(79) 



and twice that for SU(2), c.f., Eq. ([TTjl . where 72 -^ 
1 = 2/N for A^ = 2, in the hmit mg -^ 0. This dou- 
bling in the SU{2) case essentially comes about because 
for the expectation values containing terms which are at 
most quadratic in the link variables U at this order, only 
for TV = 2 we obtain contributions from two types of 
invariant integrals, Eqs. (jB3p and (|B6p . while only the 
gauge group integrations of the form in (jBSp contribute 
for A^ > 3 at this order (additional contributions similar 



to those for SU{2) here, will arise, e.g. for SU{3) at the 
next order etc.). 

For all odd TV > 3, at the present order, ^ and ^ 
therefore give. 



7n 



-( 

TV V 



1-t- 



TV^ 



JV2 



2TV 



1 

TV 



(80) 



which again leads to the same result as obtained for the 
even TV > 4 above, i.e.. 



7iv(mg) 



1 

TV 



for 



0, 



(81) 



and aU TV > 3. 

All these results are gauge-orbit independent as they 
must. While this is merely necessary, but not sufficient, 
it demonstrates that we can get away from a; = 0, at 
least perturbatively in a small x expansion. This is of 
qualitative importance as a non-zero value of a; = t/V^, 
no matter how small, corresponds to a large but finite ^ 
at i = 1 and thus eliminates the gauge freedom. 

In summary, the gauge-orbit independence of the con- 
stant in the second order term of the small x expansion 
of fh^ (x) is a direct consequence of the invariant integra- 
tions over the gauge-group elements at each site. The 
gauge-group integrations can in fact be performed at any 
order in this Taylor expansion around x = because 
the action is independent of the gauge group there (c./., 
Eq. ([55]) ). Moreover, the invariant tensor method can be 
used to demonstrate how these integrations will ensure 
that the coefficients in this Taylor expansion are indeed 
independent of {[/} at any order in x. This gauge-orbit 
independence of the mass expansion is verified explicitly 
for the constant in the second order term of (pTjl from the 
results in Appendix |B] In particular, as we have shown 
above, in D Euclidean dimensions the gauge-parameter 
expansion of the Curci-Ferrari mass in Eq. (|6ip becomes. 



(82) 



m^(x) ^ml (l - D -iN(ml) x^ + 0{x^)\ 



where 72 and 73 for SU{2) and SU{3) are explicitly given 
by Eqs. ([7T|) and fT^ . Moreover, at leading order in the 
Curci-Ferrari mass parameter, 77V (w§) is finite and of 
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the order 1/A^ in the hmit TOq 
gauge groups we found, 



lN{fhl) 



1 

TV 



0. For general SU{N) around the Gauss-Bonnet Umit of 1/^ ^ is possible 

with ^-independent mass parameter m^, to show that 
one can meaningfully define a limit TOq — > in the spirit 
of I'Hospital's rule. In this hmit, -/n -^ l/N, for SU{N) 
gauge theory with A^ > 3 (and 72 — > 1 for SU{2)). 



1 + Oimt), N^2, 
0{ml), N>?,. 



(83) 



The leading x-dependence of the Curci-Ferrari mass is 
therefore order N suppressed. Without need to adjust 
the Curci-Ferrari mass parameter fh? with x, the gauge- 
fixing partition function .^„gf of the massive Curci- 
Ferrari model therefore becomes gauge parameter inde- 
pendent at least up to the order x^ in the large A^ limit. 



VII. SUMMARY 

We have formulated the Curci-Ferrari model on the 
lattice. In the massless case this model provides an ex- 
plicit demonstration of the topological origin of the Neu- 
berger 0/0 problem of lattice BRST. The starting point 
of Neuberger's original argument was the observation of 
uncompensated Grassmann integrations producing a zero 
result in a certain limit. This same limit in the massless 
Curci-Ferrari model with its double BRST symmetry and 
quartic ghost self-interactions corresponds to the Gauss- 
Bonnet limit, ^ — > 00, of a topological model that com- 
putes the Euler characteristic of the lattice gauge group 
which vanishes for compact Lie groups. The fact that 
the Neuberger zero is independent of this special limit 
then follows directly from the gauge-parameter S, inde- 
pendence of the topological model. 

Introducing a Curci-Ferrari mass term regularizes the 
Neuberger zero. The analogue of the Gauss-Bonnet- 
Neuberger limit here corresponds to the gauge-parameter 
^— >oo limit together with the Curci-Ferrari mass rn? — > 
such that the product rn^^/S, remains finite, i.e., vr?^f% —^ 
TV? for some finite m^. In this limit, the partition func- 
tion of the massive Curci-Ferrari model on a finite lattice 
is obtained as a polynomial in the new mass parameter 
rh^ and is hence non- vanishing. The 0/0 problem is thus 
avoided. However, the massive Curci-Ferrari model is no- 
longer a purely topological model. BRST and anti-BRST 
are explicitly broken by the Curci-Ferrari mass. The 
s^(2, IR)x double BRST algebra of the massless Curci- 
Ferrari model is decontracted into a simple superalgebra 
for m^ ^0. As a result of this BRST breaking, mean- 
ing that neither BRST nor anti-BRST transformations 
are nilpotent anymore, the gauge-fixing partition func- 
tion of the massive Curci-Ferrari model is a priori not 
independent of the gauge parameter ^. This implies that 
the Curci-Ferrari mass has to depend on ^ so as to re- 
store total ^-independence, a requirement which in turn 
allows to determine this ^-dependence of m^. A gauge- 
orbit independent Taylor series expansion of m^ (^) , in D 
dimensions of the form, c./., Eq. 



.^(0^f(i 



Dji, 



(84) 



VIII. CONCLUSIONS AND OUTLOOK 

The explicit BRST breaking by the Curci-Ferrari mass 
term is well-known to lead to unitarity violations in the 
continuum quantum field theory. The BRST-cohomology 
construction of a positive physical Hilbert space breaks 
down together with this BRST breaking. There are 
explicit examples of negative norm states mixing into 
what would otherwise be defined as the physical sub- 
space [3, Q . For that reason, the parameter m^ should 
not be interpreted as a physical mass but rather only as a 
regulator for the Neuberger 0/0 problem of lattice BRST. 

In the Landau gauge limit with ^ = 0, for example, 
it has the effect of reweighting different Gribov copies 
depending on their value for the Morse potential on the 
gauge orbit, T4/[5], which avoids their perfect cancella- 
tion. The analogue of the Poincare-Hopf theorem for the 
Euler characteristic of the lattice gauge group, Eq. (|^ 
for finite m'^, becomes, 

Z„oF(t,0,m2) = (85) 

J2 sign (det Afpp ([/»*" )) exp { - mh Vu[g^'^]} • 

copies {g(')} 

At finite m^ this leads to a suppression of Gribov copies 
outside the fundamental modular region, i.e., a suppres- 
sion of all copies relative to the absolute minima of Vjj [g] . 

In particular, because we furthermore have Vulg] ex 
— trMpp, each positive (negative) eigenvalue will increase 
(suppress) the weight of a given copy. The situation is 
thus similar to that in the Gribov-Zwanziger approach 
[3|, |27J | which includes a horizon functional to suppress 
Gribov copies with negative eigenvalues, outside the first 
Gribov region. And as in the massive Curci-Ferrari 
model, this leads to a certain BRST breaking in the 
Gribov-Zwanziger framework also. While the renormal- 
izability of this framework is maintained [27], the uni- 
tarity violations when attempting a BRST cohomology 
construction of a physical Hilbert space remain to be a 
problem there as well. In fact, the horizon condition 
might be interpreted in a way as to generate a kind of 
Curci-Ferrari mass dynamically. 

The so-called soft BRST breaking of the Gribov-Zwan- 
ziger framework, which really means that it is a non-pcr- 
turbative BRST breaking, should not be a problem, if it 
is effective in restricting the Landau gauge to the funda- 
mental modular region, the set of absolute minima along 
the gauge orbits [2g| . This should be intuitively expected, 
of course, for a perfect gauge-fixing, but because of the 
existence of many degenerate absolute minima it is not 
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mathematically obvious. Note that we could at least for- 
mally achieve the same restriction here, if we define the 
Landau gauge limit of the massive Curci- Ferrari model as 
the limit of vanishing gauge-parameter, ^ — *■ 0, at finite 
mass parameter m^ = m}\f\ (analogous to the opposite 
limit ^ — > cx) at fixed m^ considered above), as this would 
imply w? — > oo in (|85p . so as to suppress all copies but 
the absolute minima of Y\j [g] . 

The arguments for the Gribov-Zwanziger framework to 
achieve this further restriction to the absolute minima in 
lattice formulations [23| , furthermore involve the thermo- 
dynamic or infinite volume limit, in which the common 
boundary of the fundamental modular region and the 
first Gribov region dominate the minimal Landau gauge 
configuration space j30i] . To find absolute minima in lat- 
tice simulations is not feasible for large lattices as this is 
a non-polynomially hard computational problem. When 
sampling the local minima of the first Gribov region on 
a finite lattice, as is usually done in minimal lattice Lan- 
dau gauge implementations numerically, BRST breaking 
and unitarity violations could therefore be a potential 
problem there also. 

In fact, related to this potential problem of minimal 
Landau gauge implementations on finite lattices might 
be the question of a gluon mass and the infrared behav- 
ior of the gluon and ghost propagators in Landau gauge. 
While early lattice studies of those propagators [3l|, [S^l 
supported the predicted infrared behavior based on their 
Dyson-Schwinger equations qualitatively well [33, |34| . 
small but significant differences are increasingly being ob- 
served nowadays in detailed comparisons and studies of 
finite- volume effects [33, |3a, 123, |3^. In particular, the 
Kugo-Ojima confinement criterion of local quantum field 
theory seems now violated in the infinite- volume limit: 

Based upon the assumption that BRST-invariance is 
unbroken, the continuum prediction is that the infrared 
dominant correlations are mediated by the Faddeev- 
Popov ghosts, whose propagator is infrared enhanced in 
Landau gauge, while the gluon propagator is found to 
be infrared suppressed. This infrared behavior is now 
completely understood in terms of confinement in QCD 
H, [3^, it is a consequence of the celebrated Kugo-Ojima 
confinement criterion. The subsequent verification of this 
infrared behavior with a variety of different functional 
methods in the continuum meant a remarkable success 
[13, SH, m. In fact, it is directly tied to the validity 
and applicability of the framework of local quantum field 
theory for non-Abelian gauge theories beyond perturba- 
tion theory. Consistent with the conditions for confine- 
ment in local quantum field theory, these predictions all 
lead to a conformal infrared behavior for gluonic Green's 
functions which is yet to be observed in lattice simula- 
tions. Because simulations must necessarily be done in 
a finite volume, where such a behavior can strictly never 
be observed, finite-volume effects have long been blamed 
for the remaining mismatch with the continuum studies. 
Recently, these effects have been analyzed carefully in 
the Dyson-Schwinger equations [33, [3a]. These results 



together with latest lattice data on huge lattices [37], |38| , 
corresponding to physical lengths of more than 20 fm in 
each direction, clearly rule out finite-volume effects as 
the sole origin for the observed discrepancies. 

Rather, the most likely origin for these discrepancies 
potentially hints at a much more profound problem: a 
BRST breaking that a sampling of minima of the gauge- 
fixing potential in lattice simulations might bring about 
much like the restriction to the first Gribov region does in 
the Gribov-Zwanziger framework. The observed evidence 
of a massive infrared behavior of the gluon propagator in 
the infinite volume limit in current lattice Landau-gauge 
implementations in fact suggests that. Any reweighting 
of Gribov copies, inside or outside the first Gribov region, 
should correspond to a BRST breaking procedure anal- 
ogous to the introduction of a Curci-Ferrari mass. Only 
if the numerical procedure converges towards a correct 
sampling of the fundamental modular region in the infi- 
nite volume and continuum limits, can the BRST break- 
ing effects be expected to go away, which they will then 
probably do together with the observed gluon mass. 

In interesting alternative procedure based on stereo- 
graphic projection to define lattice gauge fields is pro- 
vided by the modified lattice Landau gauge of Ref . [4j| . 
This gives rise to a manifestly BRST invariant lattice for- 
mulation. The Neuberger 0/0 problem is avoided there 
because it is not the vanishing Euler characteristic of 
the lattice gauge group that is calculated by the gauge- 
fixing partition function in this case, but that of a stere- 
ographically projected manifold. In this approach the 
pure lattice-artefact Gribov problem of compact C/(l) 
is avoided because the Faddeev-Popov operator is posi- 
tive, and there are thus no cancellations between Gribov 
copies |43| . Consequently, there are none along the max- 
imal Abelian subgroup [/(1)^~^ of SU{N) either. This 
is just enough, however, to remove the complete can- 
cellation of Gribov copies in SU{N) also. The remain- 
ing cancellations between copies of either sign in SU{N), 
which persist in the continuum limit, are necessarily in- 
complete because the Euler characteristic of the coset 
manifold is non-zero. It is essentially determined by that 



of the even dimensional spheres S x S 



■S' 



2N-2 



more precisely, of the corresponding even-dimensional, 
real projective spaces RP^", of one dimension less than 
the odd-dimensional spheres of the group manifold. The 
perhaps most promising feature of the modified lattice 
Landau gauge is, however, that it provides a way to per- 
form gauge-fixed Monte-Carlo simulations sampling all 
Gribov copies of either sign in the spirit of BRST. Bridg- 
ing the gap between gauge-fixed Lattice and continuum 
studies this should, in particular, resolve the present dis- 
crepancies observed in the QCD Green's functions once 
and for all, and at the same time put our theoretical 
knowledge of QCD and all gauge theories of the Standard 
Model on solid ground in a completely non-perturbative 
manner. 

Meanwhile, the decontracted double BRST supersym- 
metry of the massive Curci-Ferrari model on the lattice 
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provides an interesting testbed with a controlled BRST 
breaking and regularized Neuberger zeroes. It might have 
its own interesting topological features and interpreta- 
tion as indicated by its potential gauge-orbit and gauge- 
parameter independence. This certainly deserves further 
study, especially with regard to its potential to gener- 
alize the topological field theory relation between the 
Gauss-Bonnet and Poincare-Hopf theorems [2J, H, [4^ . 
In fact, an extension to Morse theory similar to ([55]) . 
but corresponding to the introduction of an imaginary 
Curci-Ferrari mass parameter (i.e., using a real param- 
eter (j) = im?t in (ISS]) ). and its relation to a generalized 
Gauss-Bonnet theorem reminiscent of ([5^1) were discussed 
in Ref. [46] . This might help to understand how it is pos- 
sible to achieve the gauge-parameter and orbit indepen- 
dence in the massive Curci-Ferrari model in general. And 
it might provide an interesting new topological interpre- 
tation of the Curci-Ferrari model within the decontracted 
double BRST osp(l|2) superalgebra framework. 
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APPENDIX A: BRST DERIVATION OF 

FADDEEV-POPOV OPERATOR AND 

GAUGE-FIXING ACTION 

In this appendix we show explicitly how the double 
(anti-)BRST variation in the gauge-fixing action of the 
topological model in Eq. p^ leads to the explicit form 
given in (fT4|) . This lattice transcription of a well-known 
continuum procedure is mainly given for the reader's con- 
venience and ready reference. 

Performing the anti-BRST variation on the r.h.s. in 
Eq. (fT2|) first, we obtain 



2p 



I ]~l 



F^{US 



= -J2cfFtiUa), where 

i 



is used in the standard gauge- fixing condition of covariant 
gauges. In the continuum limit it reduces to 

5^%2g^^ga^0^^4)_ (A3) 



With Eqs. (O, (flUl) we furthermore have 



(A4) 



and therefore, for the gauge-fixing action, the alternative 
form 



«C >a^ 



^s(c-(^"([/^)-h|6: 



(A5) 



As in the continuum formulation, in this form it looks 
exactly like the gauge-fixing action of standard Faddeev- 
Popov theory for the linear covariant gauge. The spe- 
cific features of the ghost/anti-ghost symmetric frame- 
work show when working out the remaining BRST vari- 
ation. From the first term we have (i), 

E (•<) ^— - 2^ E E t'^ Huf, - m (A6) 

i i j'-^i 

+ iEEtr({^.,^.}(t/S-t^))- 



4p 



I j^i 



Herein, the first term on the right implements the gauge- 
fixing condition as in standard Faddeev-Popov theory. 
The second term, containing the anticommutator {c, c}, 
is characteristic of ghost/anti-ghost symmetry because it 
combines with the remaining quadratic ghost terms to 
produce a Hermitian Faddeev-Popov operator (for any 
gauge parameter ^). To see this explicitly, consider (zi). 



T^^i^F:) - ^EEt'^(s^^^f^^- (A7) 

so that the difference {i) - (ii) yields 

EK^t^'') - -^EEt'-Nt^^-^i^)) (A8) 



2p 



I j~i 



?: E E t^ (^«f^^^ - ^^f^^:' - [^- ^4(^ + ^)) 



2p 



Y^h-F- +Y.~^1M..'^c), 



1,3 



which defines the lattice Faddeev-Popov operator Mpp of 
the ghost/anti-ghost symmetric Curci-Ferrari gauges. 
Note that the terms in (IA7I) can be written in the form 



Y.^^{sFt) = (A9) 

i 



-Ftr 



This yields the standard Faddeev-Popov operator of the 
linear covariant gauges on the lattice. It differs by the 
quadratic ghost terms in (|A6|) from the ghost/anti-ghost 
symmetric one, M-pp in (jASp . These terms are of the form 



i-5:5:tr({c-.,ca([/^.-c/p) = 



Ap 



(AlO) 



I 'J~t 



^J2^ajabcpc^^g^^b 



In lattice Landau gauge, where F^(U^) = 0, the forms 
from Eq. (|A7[) and (|A8[) both, of course, lead to the same 



14 



Faddeev-Popov operator, as they do in the continuum 
where the standard and symmetric versions differ liy an 
analogous term ex f'^'"^d^A'^/2 which vanishes for ^ = 0. 
For ^ 7^ 0, on the other hand, the two Faddeev-Popov op- 
erators do differ and, in particular, only the ghost/anti- 
ghost symmetric framework based on (jASp leads to a Her- 
mitian one which can be written in the alternative form. 



E^°^-4 



V 3 



4p^ I 



(All) 



I.J 



['^X + fL ~ '^x){'"^X + ii ^ ^x) 

+ tr([X^X''](t/,^,^-t/|,t^)) X 



\'^x\^x+ii ^x) \'^x 



-HA 



^xi^x) f 



We have added and subtracted the term proportional to 
c%c^^ in the last line here to underpin that in the contin- 
uum limit the Mpp herein reduces to the ghost/anti-ghost 
symmetric Faddeev-Popov operator. 

To complete the derivation of the gauge-fixing action in 
the ghost/anti-ghost symmetric framework, we further- 
more need work out the BRST variation of ss{c°'c°-) = 
s{c°'h°') from (O-®. This, however, is done in exactly 
the same away as in the continuum, the result is (iJi), 



5(c''6'') = 6"6" + i (c X cf 



(A12) 



Putting together all terms from (i) to {in) we obtain the 
full gauge-fixing action with extended double BRST in- 
variance on the lattice in the form, 

^OP - E { ~ *^"^"(f^') -^ctM,.nc] (A13) 

i 

+ f^"^^ + f(c. xc,f}, 
where we introduced the short-hand notation that 

3~i 

-Ktr({X«,X''}(t/^+C/j^))(c^-4)}, 

which corresponds to the ghost/anti-ghost symmetric 
Faddeev-Popov operator in (|Alip . in particular, we have 

^ -c'i Mpp^ [c] = ^ c^ Mpp^^ c] , (A14) 



^,0 



ah 



with Mppj in a simplified alternative form given by 



M, 



^5]{Retr({X^X''}C/ii)*,, (A15) 



2 Re tr 



The Faddeev-Popov operator of lattice Landau gauge was 
first derived in [13]. It might be worth recalling that 
the derivation presented there, based on the differentials 
of the gauge-fixing potential Vjj [g\ along one-parameter 
subgroups of the SU{N) gauge group, by construction 
leads directly to the Hermitian Afpp in (JAISP , and not to 
that of standard Faddeev-Popov theory on the lattice fol- 
lowing from (jA9[) . They are equivalent in Landau gauge, 
of course. Their subtle difference needs to be kept in 
mind, however, when attempting to implement smeared 
covariant gau ges for ^ 7^ on the lattice as was done, 
e.g., in Refs. [Tj, |T^. It reflects the different symme- 
try properties of standard Faddeev-Popov theory and the 
ghost/anti-ghost symmetric framework for ^ 7^ 0. 



APPENDIX B: EXPECTATION VALUES IN THE 
'WOULD-BE-GAUSS-BONNET' LIMIT 

Here we describe the explicit calculations to derive the 
results for the expectation values in Eq. ([55]) . as sum- 
marized in Table [H which are needed at the 2"^^ order in 
the expansion of the Curci- Ferrari mass parameter fh?{x) 
around a; = 0, the 'would-be-Gauss-Bonnet' limit. 

For the first term we use (fe°fe?)m2 = S'^^'Sij to obtain. 

Remember the explicit form (J15p of the gauge condition, 
FtiU^ = - ^ tr (X^igjU^.g, ~ g]u,,gd) . (B2) 

All terms in the sum of the squares of this condition 
are quadratic in g and in g^ . For A^ > 3 the only non- 
vanishing results of the group integration arise in terms 
where all the gi's match up pairwise with gj^s at iden- 
tical sites (the special case of A^ = 2 will be discussed 
separately below). We can then use for the fundamental 
g's at that site i, 



Vn^ / dgi{gi)ki {ghmn = -rjhnS, 



N 



(B3) 



In one particular term F^F"" , without summation over i, 
only the two mixed terms then survive and we have. 



V, 



N 



dg.F^F^ = -2V^^ dg, x 



(no sum i) 



Y,^v{X^glU,,g.j))(Y,^v{X'^glU, 



k'-^i 






tY{X''X''glUk^U^Jg.j) 



^-W) T.^'(9lUHU,,g,). (B4) 

The integrations over the gauge group elements at all 
sites k ^ i, attached to site i, by the same argument 
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yield a non- vanishing result only if j — k in the double 
sum over all neighbors of site i in (|B4[) . For those terms 
the group integration yields tiUjiUij — N. Therefore, in 
D dimensions, 

( {F, F) )^_^, = -^Z n '^f)^ {F, F) = (#sites) x 4D C( , 

(B5) 
for SU{N) with A^ > 3, where C^ = ^N ~ 1/N) is the 
value of the quadratic Casimir operator in the fundamen- 
tal representation. 

For A^ = 2 we obtain an additional contribution to 
{{F, F)),n2_Q from the squares of the two terms in the 
gauge condition (|B2p . This is because for g G SU{2), 



^2 ^ / dgt {gi)ki {9i)r 



■ ^krn^hi 



(B6) 



Again, however, only the squares of the same links con- 
tribute in the double sum over the neighbors of site i. 
For these, e.^.. 



(trX^g]C/,,.9,)^ 



— > - det (7,1 = Cn 



(B7) 



There are two such terms in F^F^ for each of the 2D 
links attached to site i. Therefore, the total additional 
contribution from those terms in SU(2) equals that from 
the mixed terms calculated for general N above, we have, 

((F,F))^,^= iifsitcs) X 8D Ci , for 5(7(2). (B8) 
Next, to compute the expectation value oc ( V^ )m2 o of 

Vuig] = ( E ^'-aH.gj) ( E ^^alUkigi) , (B9) 

we consider one term for fixed (neighbors) i and j in the 
double sum. Integrating this term over gj via (jB3[) . 



V^^J dgj irg\U^.jgj [^ ^ ^^dlUkigi 



(BIO) 



k, l~k 



-E ^"c glUijUjigi , 



N 



i~j 



for A'' > 3, where the only contribution arises when k = j. 
Then, I must be one of the neighbors of j. Successive 
integration over gi singles out that neighbor of j with 
/ = i, 

^N^ / ^9^J;^Y1 ^^alUijUjigi = — trC/yC/j, = 1 . 

We obtain one such contribution for every one of the 2D 
neighbors j ^ i at site i, thus summing ^^ - ^ yields 



{VS[9])fn-,o = (#sitcs) X 2Z? , iV>3 



(Bll) 



Note that the number of sites in all these expectation 
values cancels with that in ((57|) when computing the ratio 
of Eq. (P|l. 

Again, for A^ = 2 in SU{2) we obtain an additional 
contribution from (jB6p . Starting again from the contri- 
bution to (jB9p for fixed neighbors i and j as in (jB10[) we 
now obtain for N = 2, 



V^^ j dg,trg\U,,g,(^Y. ^""alUkigi) = (B12) 

2 E ^^dlUtjUjigi +^^esu trt {g\U^j)sr{9iUkj)ut ■ 



l~j 



k~j 



The second group integration over gi then produces, in 
addition to the above, an according contribution from 
the second term, which is given by 

Kj"^ / dgi -'^€su£rt {9lUij)sr(.9[Ukj)ut = (B13) 

k^j 
T ^su ^rt ^su ^vw \^ij )vr\^ij )wt ^^ T 4 CtCt U ij — i . 

This equals the first term obtained for all A^; and together 
the two again give for SU{2) twice the result for A^ > 3 
in (JB11[) above, i.e., 

( V^[g] )if,2^o = (#sites) X 4Z? , A^ = 2 . (B14) 

The hardest task here is to compute the last remaining 
term in (j63|) . ( (ic, Afppc) ) -2 n- For a start, we first note 
that 



{^c^c';^clcf),^,^^ = {S'^%, S^%i - S^^iS'^S.k) A 



ac 

ik I 



where 



Al^ ^ {{z-ccTA^-cc)]) 



3/fh^,0 ' 



(B15) 



using the notation {iccji = ic^Ci without implicit sum- 
mations over a and i, here. The expectation value in 
(|B15P is of course independent of {U}. It depends on the 
site indices only in that we need to distinguish whether 
i = j or not, i.e., 



Aab _ 



pab(ff,2\ pab _ J Pn : " 7^ ^ . 
^N V"'0) 1 ^N — 1 Q a — b ' 



l=J 



QN(fn^) , independent of {a, fe} ; i ^ j . 



(B16) 
where both P^'' and Qn are site-independent. We have 



QN{ml) 






(B17) 



from ([5S)) . while P^ has a more complicated structure, 
in general. Only its sum simplifies. 



Y.p%\fhi) ^ 2 J2 PNif^l) - J 

a, b a, fc>a 



N{ml) 
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which foUows immediately from its definition via (|B15[) 
and ((BT6)) and with Eqs. ([Mll-dSll). With these results, 

((zc,Mppc)')-,o = (B19) 

i,j>i,a,b 

+ 2 ^ <(a2) ( (Af ""Mi^" - Mt^'M^t) )rn^^ „ ■ 

i,a,b>a 

The first part with the contributions from different sites 
j ^ i would be a disaster for the intended m§ -^ 
limit: Because I'p^{fh'^)/lM{fh'^) always is of the order 
I/ttiq, we have that (5jv(™o) i^ ^^ ^^^ order 1/mg. Re- 
calling that we need to divide all terms computed here 
by the expectation value in (|57|) which is proportional 
to I'j^{fhQ)/Ii\f{fhQ), this then implies that the second 
derivative w.r.t. x of m'^{x) would contain a contribu- 
tion proportional l/friQ at x = and therefore become 
infinite in the limit ttiq —^ 0. Luckily, this contribution 
turns out to be zero for all N because of a cancellation 
between the two terms in the expectation value of this 
part in (|B19p . To see this, first consider. 



Y^Mt^-Mf^ = (B20) 

4^2 ' E E ^^ *^' {9\u^k9k) Re tr {g]U,igi) 
- ^2'EE (trC^ffetrt/j^ + trU^^trUf^ + 



kr^i Ir-jj 



trU^.tr^ +trf/|,tr[/j^ 



Because i j^ j, only the first two terms in the brackets 
contribute when integrating gi via (jB3|) . For those, the 
different group integrations over gi and gj , by the method 
now familiar, then yield. 



V, 



N 



\^ 



rbb _ r,r<f^ 






dg^dg, ^M^M^^ = 20^ 2^ 2^ -^h 4, 



a.b 



k^i l^] 



After summation over j we can replace j by k here. The 
fact that in this sum we need to restrict j ^ i does not 
matter because the non-zero contributions arise for j = fc, 
where fc is a neighbor of i and thus necessarily different 
from i. Subsequent summation over the neighbors I of k 
now, then due to the second Kronecker symbol picks up 
neighbor i of site k. We have 



^ Vj,'Jdg,dg, J2mTM^^ = 2ct E E ^'^ 
= (#sitcs)xi? (iV-1)' . (B22) 



jj¥'' 



As before there is an additional contribution from (jB6[) 
for SU{2) here also. As i and j are different, after in- 
tegrating gi this time, this contribution is nonzero only 
for precisely the other two terms in the brackets above. 
Using (|B6[) it is straightforward to verify that the result 
for those two terms again matches the contribution just 
calculated from (|B3[) . We can therefore summarize that 
for all N including TV = 2, 



^ (M-Mj*) _ , g - (#sites) X 2D 2C/' (l + 6^.2) ■ 

i,j^i,a,b 

(B23) 
For the second term in the first part of (|B19p with i ^ j, 
which allowed us to drop the diagonal terms in both 
Faddeev-Popov operators, we need that in the fundamen- 
tal representation. 



{X^X"-)ij{X''X'')ki = i[N ~ j^) SaSjk + -^jp 5ij5ki , 



and 



{X°'X'')^j{X''X%i = -—-5a5]k + T (1 + TTj) 6ij5ki 



2N 



m) 



(B21) so that we can write (for i ^ j), 



J2 M^'M^^ = 4^ ^ Re tr {X''X^glu,kgk) Re tr {X-X''g}Ujigi) S^jSh 

a,b A;~?' Ir^j 



(B24) 



= EE (tr^'^^C/f, + trX'^X^(7j',)(trX'^X''C//, + tr X'X''^^ SkjSu = 
EE{^(^'-2) + ^tr[/i^tr[;--^(tr[/i^C/-+trt;-t;-)+i(l + ^)((trC/^)%(trf/^)^)}MH 



fc~Z Ir^j 



The group integrations via (|B3p and, in addition, for the special case of A^ = 2 via (jB6p proceed as before. The 
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explicit results of the corresponding calculations above These are the basic terms that arise at the quadratic 



can all be reduced to essentially using two relations sum- 
marized as follows, 

{tr {glUtjgj)tr {glUkigi))~2^^ = SuSh + SN.i^tkSji , 
i^r {glUijgjglUkigi))^2^Q = N SuSki - SNa^ik^ji ■ 

(B25) 



order in x, and hence in the gauge-transformed links U^ , 
of our mass expansion (recall that the terms linear in x 
vanish in this expansion because the linear order terms in 
U^ do upon the group integrations). Using the relations 
(lB25l) in ((B24| . we obtain, 



a, 6 



N^-2) ■ ^ 



and hence, together with (jB23|) . 



27V2 



JN,2 



ii + li 



1 + 



k^i l~j 



m 



i,j^i,a,b 



{{M^'M^~Mf^Mljt))-,„ = (#sitc..)2D X i 






(B26) 



(B27) 



The two terms from (jB23p and (|B26[) therefore cancel 
and the first part in (|B19p thus vanishes in either case, 
whether A^ = 2 or A^ > 3. At the same time this can- 
cels the otherwise quite disastrous singularity of the mass 
expansion in the frig -^ limit, as promised. 

In the second term in (|B19p . the one with i — j, we 
need products of diagonal entries of the Faddeev-Popov 
operator of the form (c./., Eq. (jA15P : no implicit sum 
over i here either), 

^^' = - E ^ i^' {^"' ^'i ^'k + tr {^^ X"} f/,y 



k.r^i 



E{5^^"'(trf/4 + trf/|, 



(B28) 



+ \d'^'^{ir^X^Uf, + ir^X'^Ul)] 



Where we have used the identity 



{X^^X^} = -^5"^ -idf^^^X" 



(B29) 



For SU{2) we set d"''^ = which then leaves us with only 
the first term in (|B28|) . Because, e.g., 



(tr {glU^kgk) tr (iX'^gjUug^)) 



(B30) 
= 0, 



due to the tracelessness of the generators AT", there will 
be no terms linear in the AT's in the expectation values 
of squares of the operators in (|B28P for A^ > 3 either. 
The terms quadratic in the X's simplify from relations 



as follows, 
{tv{iX-glU,kgk)tr{iX''g}Uug,))-^2^^ = (B31) 

^{tr {tX^gjUuU.kgkiX'^)), 



N 



' m2 



^ S^'Ski. 



2N 



We do not have to worry about corresponding terms from 
(IB6|) for A^ = 2 here because these contain the d-symbols 
which are zero in SU{2). With (|B25)) . (|B30)) and (|B3l1) 
we therefore obtain, 

(Mf,^M^'^),_^,^= (B32) 

We therefore have, implicitly summing over color indices 
again but not over sites i (yet), 

( iic'^M^''c% (ic^M^ V),) _ ., „ = (B33) 



N^2 2D_ ( 2 xabxcd 

4N 



2D 



(r'^ic^c^r'^c-'c^)^,^^ 



2D 

w 



((^xc)?)rn=,o 



for A^ > 3, and twice that for N — 2 again, where the 
d's are zero and where the x = expectation value of 
the remaining first term above agrees with that of the 
quartic ghost interaction. 



lUff^l) 



((^-x^)?)„..o=((-"^")?)^^.o = i^ 



(B34) 
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Finally, summing in (|B33|) over the sites i of the lattice, 
and from (|Bl9)) with ((B26)) . we obtain 



((ic,Mppc)')-, = 2I?(#sites)x 



AN 



((c X c) 



;»Ti2 



_2_ Jiv(mg) ^ 



N iN(ml) ' 
1 /^'(m§) 



(B35) 
>3, 



'4<(""-")'>"^-°-4/.(mg)' 



7V = 2. 



The only piece left to compute is the expectation value 
(at a: = 0) of the quartic ghost interaction, 



_ JN{ml) 






(B36) 



The integral expression for Jn , analogous to that for In , 
c./., Eq. ([38)l . is given by 



.N-l 



JN{m^) = / n d{ic"')dc'' (-i(cxc)2) (B37) 

•' a=l 

exp < JTO^ c-c — o (c X c)^ > . 



For SU{2) and S't/'(3), respectively, the resulting J2{rn?) 
and J3(m^) are given explicitly in Eqs. ([M]) and ((57)) . 
This completes the computations of the terms in (|63p . 
The results are summarized in Table HI 
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